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EVALUATION OF THE PROBABILITY INTEGRAL 

Br Fbank Oilman 

A simple and convenient formula for calculating the numerical value of 
the definite integral 



/' 



-**d« 



has long seemed to me a desideratum. The three series commonly used for 
computing this integral are not convenient in practice, except when the argu- 
ment x is either very small or very large. It is believed that the formulas 
here given will be found equally convenient for any value of x from to 3, 
and that they will, in general, give the integral correct to five figures. It is 
true that some of the old tables give this integral to eleven, and some even to 
fourteen decimals ; but so great a degree of accuracy appears to be unnecessary. 
Besides, the method here given can be extended so as to secure any degree of 
accuracy desired. 
The integral 



X 



e-*'dx 
'o 



is of great importance in many branches of applied mathematics. It is used 
in discussing the effect of refraction,* in investigating the secular cooling of 
the earth f and in discussions in regard to the conduction of heat, as found in 
Fourier's work, Thiorie anaiytique de la chaleur. 

*See Chaurenet's Astronomy, vol. 1. 

t See Thompson and Tait's Statural Philosophy, vol. 1, p. 717. 
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It is of the greatest importance in the theory of probability, or proba- 
bility of eiTors, and on this account, when written in the form 



r 



e - ** dx = erf a; 

'x 



it is sometimes called the error function of x, while 

er*'dx — erfcx 



is called the error function complement of x. The former of these is des- 
ignated by the notation erf x, and the latter by erfc x. The two are so related 
that 

erf x + erfc x = / tr* dx, 

or 

erf x + erfcse = JyV. 

We shall give formulas for computing directly erfc x, since, this being 
known, erf a; can be immediately deduced. 

Besides the applications of erf a; above mentioned, it is of great use in 
the integration of many differential equations (Kramp says an infinite number) 
and Glaisher in the Philosophical Magazine for 1871, * has given a list of 
some of the definite integrals which can be expressed in terms of it. 

Glaisher says that from its uses in physics, erf x may fairly claim at 
present to rank in importance next to the trigonometrical and logarithmic 
functions. 

The following is the formula proposed for computing erfc x for values 
of x from to 1 : 



i: 



e~ x ' dx = A + ax + bx*+ ex 6 . (1) 



When x > 1 we may write 



I 



"e-*'dx = A + -+ ^ + - B . (2) 

X X s X 6 v ' 



•Vol. 42, pp. 294, 421. 



1908] 



EVALUATION OF THE PROBABILITY INTEGRAL 



The following are the values of the constants A, a, b, and 
ing to assigned limits of x. 



169 
c, correspond- 



Limits 


of X 


A 


Log a 


Log b 


Log c 


to 0.40 








9.5221346w 


8.9674202 


0.40 ' 


' .55 


0.000225 


9.9995025 


9.5148674m 


8.9090072 


.55 ' 


« .70 


.002127 


9.9966907 


9. 4956431m 


8.8260698 


.70 < 


' .85 


.009435 


9.9882445 


9.4598519ft 


8.7203776 


.85 ' 


' 1.00 


.024000 


9.9739908 


9. 4159907n 


8.6187244 


1.00 ' 


' 1.15 


.960956 


8.8864253m 


9. 2954356m 


8.7802668 


1.15 « 


' 1.30 


.871946 


9.0553760 


9.5600979m 


9.0974426 


1.30 ' 


< 1.60 


.809877 


9.4213859 


9.7244149m 


9.3195134 


1.60 ' 


' 2.00 


.831292 


9.2883208 


9.5942662m 


8.9499540 


2.00 ' 


« 2.40 


.883119 


7.5394525 


8.9915539 


9.6821021m 


2.40 < 


' 3.00 


.893330 


8. 6048101m 


9.3971250 


9.8547853m 



The following was the method of calculating these quantities. Assume 

a' + b'x* + c'x* = e~ x \ 

Substitute in this equation certain values of x, and from the condition 
equations thus obtained, form the normal equations, according to the method 
of least squares. The solution of these equations will give the values of a', 
b', and c'. Then by integration we obtain 



i: 



e~ x, dx — A + a'x + -^ x 3 + 



in which A is a constant, a — a', b = J6', c = \c'. 
When x > 1 we write 

a' V c' 

X* + X* ^ X 6 ~ 

and proceed in the same way, remembering that a = — a', b = — J6', c = — .£c'. 
In order to explain more fully the method in detail, let us find the coeffi- 
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cients A, a, 6, and c corresponding to the limits of x : 1.60 to 2.00. Substi- 
tute in the equation 



^ b - *L 

X* X* X 6 



-ar> 



the following successive values of x: 1.6, 1.65, 1.7, 1.75, 1.8, 1.85, 1.9, 
1.95, 2.00 and write the 9 equations of condition : 



0.39062496a' + 
•36730950o' + 
.34602080o' + 
.32653068o' + 
.30864200a' + 
.29218412a' + 
.27700833a' + 
.26298491a' + 
,25000000a' + 



0.152587866' + 
.134916306' + 
.11973039*' + 
.106622296' + 
.095259876' + 
.085371556' + 
.076733616' + 
.069161056' + 
.062500006' + 



0.05960463c' = 
,04955605c' 
.04142920c' : 
.03481545c' : 
.02940120c' : 
.02494421c' : 
,02125585c' = 
,01818831c' = 
,01562500c' : 



0.07730473, 
.06571026, 
.05557621, 
.04677061, 
.03916391, 
.03263074, 
.02705186, 
.02231492, 
.01831564, 



where the second member is in each case the value of e -x *. 

From these condition equations the normal equations are formed in the 
usual way, and are as follows : 

0.90288292a' + 0.294819906' + 0.09812945c' = 0.12839856 
0.29481990a' + 0.098129456' + 0.03324910c' = 0.04358243 
0.09812945a' + 0.033249106' + 0.01145041c' = 0.01502736 

Solving these equations we find logo' = 9.2883208n, log 6' = 0.0713875, 
log c' = 9.6489240ra ; and from the known relations between o', 6', c' and 
a, 6, c, we have logo = 9.2883208, log 6 = 9. 5942662m, logc = 8.9499540, 
as given in the table. The constant A is determined from the condition that 
when x = 1.6 formula (2) should give the correct value of 

/1.6 
e-* ! dx. 

A similar method was used in obtaining the other coefficients of the table ; 
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but the number of condition equations was not always the same, and in some 
cases only 4 were used. 

Theoretically it would seem as if a large increase in. the number of condi- 
tion equations would add largely to the accuracy of the results, but in practice 
I have not always found this to be the case. 

In order to illustrate the use of the table let it be required to find the 
value of 



/•1.61 

/ e 



e~ x2 dx. 



The computation is as follows : 

loga = 9.2883208 logft = 9.5942662n logc = 8.9499540 

log 1.61 = 0.2068259 log (1.61)» = 0.6204777 log (1. HI) 6 = 1.0341295 

7.9158245 



9.0814949 


8.9737885n 


A = 


0.831292 


ax = 


0.120641 


ex* = 


0.008238 




.960171 


bx 3 = - 


- .094143 



.866028 



*dx 



loga= 9.9995025 

log .45 = 9.6532125 

9.6527150 



logi = 9.5148674n 
log(.45) 3 = 8.9596375 
8.4745049w 



logc = 8.9090072 

log(.45) 5 = 8.2660625 

7.1750697 



A = 0.000225 
ax = 0.449485 
ex 6 = 0.001496 
0.451206 
&b 8 = 0.029820 



.421386 



/.4S 
6 



>dx 



.866028 
.421386 
.444642 



/•1.61 



dx. 
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If we attempt to compute 



fi.ei 
Jo 
by the three series commonly used, namely, 



jf 



X 3 1 X s 1 X 1 

X 1- — — 4- 



3 1.2 5 1.2.3 7 r 1.2.3.4 9 

, l\ 2x * (2a: 2 ) 2 (2ce 2 ) 3 (2x 2 ) 4 -i 



* V,r 2x L 2x 2 + (2<c 2 ) 2 (2x 2 ) 3 + 



1.3.5 1.3.5.7 



(2a: 2 ) 2 (2x 2 ) 3 ' (2x 2 ) 4 J' 

we shall find that it will require 9 terms in each of the first two series to give 
the integral correct to three figures, while it cannot be obtained correct to three 
figures by using the third series, no matter how many terms are taken. 

As previously stated, it is believed that the formulas which have been 
given will, in general, give the integral correct to five figures, for the reason 
that I have made comparisons of integrals computed by these formulas, for 
values of x from to 2, and at intervals of every 0.02, with the values as 
given in standard tables. In the hundred comparisons thus made I found no 
discrepancy greater than one in the fifth figure. 

The reason why the method above explained should be expected to give 
accurate results was explained in an article on the ballistic problem, published 
in the Annals of Mathematics for April, 1905. 

It was not deemed advisable to extend the table of coefficients for values 
of x greater than 3, as the third series commonly used is very convenient for 
values of x beyond that limit. 

Boston, Mass. 



